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SUMMARY 


This paper deals with the application of system identification methods 
in flutter testing of aeroelastic structures* The accuracy with which flutter 
parameters are estimated depends upon the test plan and on the algorithms used 
to reduce the data* The techniques for selecting the kinds and optimal posi- 
tions of inputs and instrumentation, under typical test constraints, are pre- 
sented. Identification results for both the input/output transfer function 
and the values of physical parameters are presented. Numerical results on the 
optimal input spectrum and the accelerometer location for estimating flutter 
parameters of a two dimensional wing are obtained using these algorithms. 
Current work on applying system identification methods to high order three 
dimensional aeroelastle structures is discussed. 


INTRODUCTION 


The objective of flutter analysis is to quantify the critical points or 
boundaries of flutter and the stability margins associated with subcritical 
responses. While it is true that analytical predictive techniques have become 
increasingly useful to this objective, actual testing and data analysis is 
always required for verification of these analyses, or to provide results 
where analytical assumptions are suspect. Thus, flutter test analysis tech- 
niques are being developed which use experimental data (usually noisy) to 
provide accurate estimates of both subcritical stability margins as well as 
aid extrapolation to the critical points (refs. 1 and 2). To be most useful, 
these techniques should provide real time (or near real time) estimates to keep 
test times at a minimum. 

Further requirements on these test analysis techniques are emerging due 
to new aircraft concepts. New structural concepts, such as light weight com- 
posites technology, and control concepts, such as the active control of maneu- 
ver loads and flutter margins, will require multivariable testing analysis 
methods. These multivariable analysis techniques are necessary to define the 
modal frequencies and damping of many interactive structural components in 
complex aerodynamic regimes. 

To meet the challenging requirements of estimating accurate subcritical 
flutter test parameters and to use these results to effectively predict flutter 
boundaries for multivariable systems, a systematic approach must be adopted. 
This approach should integrate the specification of test instrumentation and 
inputs with multiinput/multioutput data analysis procedures. 
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The key elements of such an aeroelastic integrated testing analysis of 
a model or of a prototype vehicle are shown in figure 1* Firsts the test 
objective must be quantified* Historically^ this test objective has progressed 
from finding the flutter boundary to more current determination of the fre- 
quency and damping of the subcritical stability margin. The need to be able to 
better use subcritical data to predict the boundary requires determination of 
the parameters of a flutter model which may contain two or more states of the 
system* Of course, accuracy specifications for these various levels must be 
set. Second, the operating points (of a wind tunnel or flight regime) must be 
set to provide the basis for meeting the objectives within test safety con- 
straints. 

To implement the test objectives at the required points, an extensive 
analysis of test inputs and instrumentation will minimize the probability of 
ineffective results due to the improper excitation of critical modes and low 
signal/noise ratios. With the test configuration specified, the data are 
collected and analyzed using either a spectral analysis technique (e.g., fast 
Fourier transform (FFT, ref. 3) or Randomdec (ref. 4)) or an advanced parameter 
identification algorithm. 

This paper focuses on the specification of test inputs and instrumentation. 
Specifically, the three major elements of the test configuration are: 

(a) Choice and location of instruments (e.g., accelerometers, strain 
gages, gyros). 

(b) Choice of inputs with respect to type (e.g,, sinusoidal, swept sines, 
random), and location of inputs and frequencies, and energy of 
inputs . 

(c) "Required capability of test analysis procedures. 

Analytical methods for input design and instrtmient selection to obtain the most 
accurate estimates of parameters in models describing the flutter behavior of 
aerodynamic structures are developed. The methods, based on system identifica- 
tion technology, minimize the expected covariance of errors in estimates of 
unknown parameters. The locations of the instruments and the inputs (if vari- 
able) may also be optimally selected. 

This paper describes a simple model of an aeroelastic wing. The dynamics 
of the wing can be formulated as a state variable model. The analytical formu- 
lation of the input design problem for state variable models with unknown para- 
meters is given, along with a description of the methods used for selecting 
the kind, accuracy, and locations of instruments. Some results on the selec- 
tions of instruments and inputs to accurately identify the flutter character- 
istics of a two dimensional wing are described. Finally, the techniques are 
applied to large aerodynamic structures and the conclusions drawn from this 
work are discussed. 
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STATE SPACE EQUATIONS FOR A TWO DIMENSIONAL WING IN FLUTTER 


Flutter is an interaction between nonsteady aerodynamic forces and elastic 
forces in a structural component* To study the experimental design techniques 
for flutter testing^ a reference model for the flutter of a two dimensional 
wing is given based on the work of Houbolt (ref* 5)* The symbols used follow 
those of reference 5« 

An oscillating two dimensional airfoil in an incompressible flow can be 
modeled as shown in figure 2* Various forces acting on the airfoil ares 
(a) lift L^ at quarter chord and lift L 2 at three-quarter chords (b) restoring 

force and moment through the elastic axis located at (a) ^ (c) force and moment 
associated with the inertia of the substance constituting the medium (these 
will be neglected) j and (d) external forces and/or moments ^ used to excite 
flutter or inadvertently transmitted through the structure. The lifts L^ and 

L^ are modeled with appropriate delays, Houbolt (ref, 5) shows that the aero- 

elastic equations for the wing can be written in terms of nondimens ional vari- 
ables as follows (see also fig. 2)2 



By defining 

0 

4 w 
(j)^ A ^ 


145 



equation (1) becomes 
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where I is an identity matrix* This is the state space representation of an 
aeroelastic wing and can be written compactly as 

X = Fx + Gu (4) 

where x is the state vector, u is the input vector, and F and G are transition 
matrices which contain unknown parameters* An accelerometer placed at e will 
measure ^ 



w. 



(5) 


« « 

The quantities w- and (j) can be expressed in terms of x using equation (3) . 
Then ^ ^ 


y = Hx + Du 

The transfer function between y and u is 


( 6 ) 


^ - H(.I-F)-1g + D 

4 3 2 
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This transfer function can be written again into a state space form, often 
referred to as a canonical form. 


OPTIMAL SELECTION AND LOCATION OF INPUTS AND INSTRUMENTS 


As shown above, the flutter equations of a wing can be written in either 
the state variable form or the transfer function form. The multivariable state 
equations and measurement equations are equations (4) and (6). In practice, 
the measurements, y, are corrupted by additive noise, v, so that 


y = Hx + Du + V 


( 8 ) 
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where v is assumed to be a white noise source with power spectral density 
matrix R. The unknown parameters and some control parameters (locations of 
inputs and instruments) are imbedded in the matrices F, and H. The unknown 
parameters 5 whose estimated we are interested in^ will be denoted by 9. 

The accuracy of the parameter estimate 0 is expressed in terms of the bias 
and covariance properties of the estimate* It is assumed that an unbiased and 
efficient estimation procedure is used so that the input design and instrument 
selection can be carried out independently of the estimation procedure* This 
makes it possible to compute errors in the parameter estimates based on the 
Cramer-Rao lower bound* This bound is computed around an a priori value 0^ 

for the parameters 9, The information matrix M is related to the error in 
estimated by the following relation 

cov(0 - 9) > M ^ (9) 


where 9 is the estimate of 9* 

The information matrix depends upon the input energy distribution and its 
location and instrument accuracies and their locations. The design procedures 
presented here will work with the properties of the information matrix* For 
physical reasons^ a quadratic constraint is placed on the inputs and the state 
variables 



(x*^Ax + u*^u) dt 


< E 


( 10 ) 


where A is a symmetric positive semidefinite matrix. An equation for the infor- 
mation matrix^ under the constraint of equation (10) ^ in the frequency domain 
is now obtained* 


Information Matrix in the Frequency Domain 
The relation between y and u in frequency domain is 

y(w) = - F)"^G + D} u(w) 

A T(tja,e) u((jj) 


( 11 ) 
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where TCa)^©) is the transfer function from the inputs u(o))j to the measurement 

y (w) » 


Equation (10) is written in the frequency domain as 



u(oo) u*(a)) + Tr(Ax(u)) x*(o)))} do) < E 


o 


( 12 ) 


where Tr is the trace operator and denotes conjugate transpose. If S(a),6) 
is the transfer function between x and u, equation (12) may be written as 


or 



00 

{1 + Tr(A s(co,0) s*(w,0))} u(w) u*(co) dw <_ E 


o 


I 


00 

X(o),0) u(03) u*(o)) dw = E 


o 


(13) 


The inequality sign can be removed for linear systems because increasing the 
input amplitude will increase the accuracy of all parameters. The information 
matrix for parameters 0 from measurements y, per unit time, is as follows (see 
refs. 6 and 7 for details): 


M = 



R ^ u((jo) u*(w) dw 


(14) 


Defining 


u(oo) A \®(w,0) u(o)) 


(15) 
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Equations (13) and (14) become 


/ 


00 

u((jO) U*(0)) 


o 


d03 


E 


M = 



-1 ^ 1 

^ 39 X(a),0) 


u(o)) 


U* (O) ) d(l3 


(16) 


The information matrix^ M, serves as the basic quantity upon which the 
input and instrumentation requirements are to be detennined. Maximizing M by 
appropriate input and instrumentation design parameters leads to output data 
which have a high information content on the system parameters. That is, the 
sensitivity of the outputs to parameters, for example, is maximized by exciting 
the modes which are most affected by the parameters. Basing the design on M, 
though mathematically simpler, has some disadvantages in practice. If the 
trace of M (e.g., the sum of diagonal elements) is maximized, an almost singu- 
lar information matrix may result. The inverse of M is the lower bound on the 
parameter covariance matrix. If M is nearly singular, its inverse may contain 
large diagonal elements, leading to large errors in the estimates. 

For this reason, it is more desirable to work directly with the inverse 

of the information matrix, M This matrix can be viewed as the ellipsoid 
of uncertainty of the parameters. Though mathematically more difficult to 
minimize, this matrix gives useful results since we are minimizing the para- 
meter covariances directly. Two types of methods can be used to minimize 

M These are based on the following functions of M 

(1) Minimize Det (M ^): This method will minimize the volume of the 

uncertainty ellipsoid. This also minimizes maximum error in the 
estimate of the transfer function. 

(2) Minimize Tr (WM i This method minimizes a weighted sum of the 
parameter estimate covariances (W is the weighting matrix which 
penalizes certain estimate errors more heavily than others). The 
weighting matrix serves two purposes. Since the covariances of 
different parameters have different units, the weighting matrix con- 
verts each term in the sum to the same units. Secondly, the weight- 
ing matrix offers tremendous flexibility because it is possible to 
assign varying importance to parameters through weights on their 
nondimensional covariance. This is considered to be one of the most 
suitable performance criteria, since it works with parameter estimate 
covariances directly. 
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Choice and Location of Optimal Input 


The optimal input possesses certain properties which are quite important. 
They are presented here without proof (ref, 6) s 


(1) The optimal input has a discrete spectrum (or point spectrum). The 
number of frequencies with nonzero power does not exceed 

where m is the mnnber of parameters. . 

'■ "I- 


(2) If the spectrum of u contains fewer than m/ (2p) frequencies » the 
information matrix is singular (i.e., all the parameters cannot be 
identified) . 

1 ^ 

(3) The optimal input which minimizes Det(M ) satisfies a minimum output 

error criterion. In other words, this input gives the best estimate 
of the transfer function. 


(4) They satisfy two important theorems (see refs. 6 to 8), which convert 
this complex nonlinear problem into a computation technique. 

It has been demonstrated that the computation procedure summarized in 
appendix A can be applied to select the input spectrum which gives the desired 

minimum of M 


Practical considerations in the computation of optimal input .- The algo- 
rithm of appendix A will produce an optimal input design with a sufficient num- 
ber of iterations. However, at each iteration, the procedure adds one point ' 
to the spectrum of the input. For practical implementation, it is desirable 
to have as few frequencies in the optimal input as possible. During the compu- 
tation, a few steps can be taken to reduce the number of points in the spectrum. 
Suppose the normalized input at any stage has k frequencies co. with power a. 
(i=l,2,...,k). Then: ^ ^ ^ 

(a) Frequencies less than Ao) apart can be limiped into one frequency. 

* 

Suppose q frequencies are within a band Aw wide. Then they can 
be replaced by one frequency w* with power a* where 


q 

a* = S 
i=l 


a. 


and 

w* 


_ 1 _ 

a* 


q 

E 

1=1 


•k k 
a. 03. 

1 1 
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(b) From this new input, all frequencies O)^ with power less than a 

threshold are dropped. The remaining frequencies do not satisfy 
the constraint of equation (13), so the design is renormalized. 

Steps (a) and (b) should be carried out to ensure that the design does 
not become degenerate. This "practicalization" requires judgment of Aw and a. 

Choice of location of input .- The transfer functions T(w,9) (the input- 
to-output transfer function) and S(w,9) (the input- to-state transfer function) 
are both linear functions of the control distribution matrix G. The location 
of the input affects G in a linear fashion. Therefore, if 3 is an input loca- 
tion parameter, the transfer functions T(w,9) and S(w,9) can be written as 


T(w, 9) = T^(w,9) + 3 T 2 (u, 9 ) 

S(w,9) = S^(w,9) + 882 ( 0 ), 0) , 0 < 8 1 1 


(18) 


Equations (13) and (14) can, therefore, be written as 



00 

{ 1 + Tr (A (S^ (w , 9 ) + 882 (o) , 9 ) ) 


o 


(S^ (oOj 9) + 3^2 (iA), 0) ) *) } u(u)) u*(u}) du) = E 




and 


Y(1 + C ^8 + = E 

M = + 28M^2 + 


(19) 


( 20 ) 


Y is a scalar which adjusts the energy in the input to satisfy the quadratic 
constraint on the input and the states. Equations (19) and (20) can be com- 
bined into one equation. 


M = 


l-hC^3 + C23 


2 +28M^2 + 


-M22] , 0 < 8 < 1 


( 21 ) 
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B can be selected to minimize |m or Tr(WI In fact 5 an algorithm similar 

to that of appendix A for input frequency and power selection can be developed. 

In a more general case^ when there is more than one input and each input 
may be placed over any point in two or more dimensions ^ the number of para- 
meters B which must be selected optimally is more than one. The optimization 
becomes somewhat more difficulty but the basic approach remains the same. 


Choice and Location of the Instruments 

In addition to selecting the location and type of the excitation signal, 
there are two other design consideratiGns in planning a flutter test. These 
are the determination of the kind of instruments which must be used to record 
flutter response and the choice of instrument location (if there is a choice). 
Though the problem of instrument selection and location can be treated simul- 
taneously, for sake of simplicity we treat them separately. 

Selection of instruments .- The selection of instruments is a tradeoff 
between dynamic range, accuracy, and cost. The dynamic loads are often 
limited by structural constraints, and it will be assumed that the instruments 
cover this range. The accuracy with which the parameters may be estimated is 
then determined by the accuracy of the instruments. It is clear from equation 
(14) that the information matrix has an inverse relationship with the measure- 
ment noise covariances. 

00 

R ^ u(w) u*(w) dw ( 22 ) 

o 


For the purpose of instrument selection in general, the measurement noise 
covariance matrix is diagonal, i,e,, 

R"^ = diaglr^^, (23) 

where is the covariance of random noise in the ith instrument. The total 

cost of the p Instruments is a sum of the cost of individual instruments 


G = 


P 

Z 

i=l 


C.(r..) 

X XX 


(24) 


The total cost of the instrument package is assumed to be fixed. Either of the 
criteria of equation (17) may be minimized under the cost constraint and 
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X 


(25) 


> 0 
XI — ' 

The Lagrange multiplier approach may be used for optimization. For examplej 

If the criterion requires the minimization of |m ^|, the modified cost function 
is 


J = |m ^1 + + UC (26) 


where X^, i=l,2,,o,,p, and y are Lagrange multipliers. The following optimiza- 
tion equations result: 

’'ii = ° 

riC 

^ I 1 2 -1 8M . i 

* Imp Tr {M ■" - u 9 ^ = 0 , i=l,2,...,p . (27) 

ii ii 

Equations (24) and (27) are 2p+l equations in 2pH-l unknowns X^, r^^ and y. 

Note that if any r^^ is zero, the corresponding instrument has infinite error; 
in other words, this instriament should not be used. 

The optimal value of r^^ would act as a guideline in selecting the instru- 
ment. Often, it is not possible to obtain an instrument with mean square error 
of — ^ exactly and cost C.(r,.). 

r . . X XX 

XX 

^ Location of instruments . - The transfer function T(a),6) is a linear func- 
tion of the measurement distribution matrix H, and, therefore, the position of 
the instrument. For this reason, optimal choice of instrument location can 
be determined in the same was as the optimal positioning of inputs. 


RESULTS 


To demonstrate the application of the methods described above to multi- 
variable flutter problems, a two dimensional wing is considered. The values 

2 2 

of the parameters are as follows: y = 10, (k^/c ) = a = 0.1, a/c - 0.35, 

r^^ = 0,1, - 0.4, a^ = 0.6, and b^ = 0.3. The velocity is taken as 15.25c 

'^meters/sec (50c ft/sec) and the natural frequencies of the rotational and ver- 
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tical motions are 10 hertz and 2 hertz ^ respectively « This gives the system 
matrices for the nondimensional state equations as 


and 


F = 


0 

0.05 

0.1 

-0.01579 

0 

0 

-0.2 

0.1 

0 

-0.3948 

•0.6 

0.684 

-0.372 

0.01895 

0.1105 

1 

0 

0 

0 

0 

0 

1 

0 

0 

0 


(27) 


e e 

= [0.1 , , -0.12 -0.48 , 0 , 0] 


(28) 


The measurement distribution matrices are 


e e 

H = [0 , 0.5 -0.2 —, O.l + O.l-^, 

c c 


-0.01579 , -0.3948—] (29) 


A a 

D = (0.1 + ) 

e c 
c 


where 


and 


are parameters which define the locations of the input actuator 


o 

c c 

and the accelerometer. The noise in the accelerometer is assumed to be white, 
with a standard deviation of 0.02 in dimensionless tmits (this corresponds to 
about 0.61 meters/sec/sec (2 ft/sec/sec)) and a sampling interval of 4 milli- 
seconds. It is assmed that we are interested in estimating the parameters 
"2 -2 


0 ) 


r’ “(j)» ^2’ ^2’ ^2’ 


a. 


The poles and zeros of the transfer function between the measurement and 

ef e^ 

the input — and — both equal to 0.1, are, in radians per second (the nondimen- 
sional values are multiplied by 100) 
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Poles 


-5.76 ± 55. 6j 
-19.4 + 16. 6j 
-6.83 

Zeros -0.163 + 76. 2 j 
-30.0 
-5.68 


(30) 


Input Design 

As mentioned before, the inputs are designed to maximize the identif iabil- 
ity of flutter parameters from output. The test duration was selected as 2 
seconds. This is a fairly short test, but in terms of the natural frequency of 
the wing in flutter, it is long enough so that the steady state input design can 
be applied. In the design procedure, the locations of the excitation and the 
accelerometer are kept fixed, and the input frequency spectrum and power in each 
frequency of the spectrum are selected. 

To simplify the procedure, the parameters in the F, G, and H matrices are 
identified directly (D does not contain any unknown parameter). The input is 
designed to produce the most accurate estimates of the xmderlined parameters in 
equations (27) and (28) in the sense of minimizing the determinant of the in- 
verse of the information matrix. Starting from the topmost spectrtun of figure 
3 with available power distributed equally among five frequencies at 2, 4, 6, 8, 
and 10 cycles per second, the iterative design procedure gives the results shown 
in the figure. In every iteration, the design procedure adds a new frequency 
or increases the power at a frequency already included in the spectrum. This 
leads to a large number of frequencies in the computed spectrum. As mentioned 
above, this design can be simplified. When frequencies with relative power less 
than 5 per cent or closer than 0.4 hertz are merged with the neighboring fre- 
quencies, the resulting design is shown in figure 4, There are eight frequen- 
cies in the optimal spectriim — three each clustered around the two oscillatory 
modes and one each at a low and an intermediate frequency. Of the three fre- 
quencies around the oscillatory modes, one is below, one is above, and one is 
close to the natural frequency. This characteristic seems to be quite general 
and substantiates Gerlach’s intuitive choice of input frequencies for the iden- 
tification of the short period parameters of an aircraft (ref. 9). A 2 second 
time trace for this input spectrum with initial phases selected at random is 
shown in figure 5. 


Choice of Accelerometer Location 

In the flutter analysis of a two dimensional wing, typically only one 
accelerometer is used. It is generally desired that the best accelerometer 
within the test budget be selected. There may be a possibility of using two 
poorer quality accelerometers. This tradeoff was not studied. 
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The location of this accelerometer is an important parameter. The follow- 
ing performance index is considered^ 



where 6 ^ (i=l, 2 , . . . ,m) Is the set of parameters of interest and 00 ^ is the 

standard deviation of estimation error in the ith parameter. Figure 6 shows 
the values of the performance index for constant input rms value and constant 
output rms value as a function of the accelerometer position. The minimum of 
the curves gives the positions of the optimal location of the accelerometer 
under the two constraints. 


Simulation and Maximxmi Likelihood Identification 
The flutter equations for the transition matrices of equations (27) to (29) 


Both — - and — are taken as 0.1, and 
c c ’ 


are simulated with the input of figure 5 

the accelerometer rms random error is 0.02 in the nondimens ional units 
transfer function between the input and the output is as follows: 


The 


-0.0155 


(s^ + 0.360s^ + 0.599s^ + 0.208s + 0.0099) 

(s^ + 0 . 5 12s* + 0 . 45 7s^ + 0 . 158s^ + 0 . 0292s + 0 . 00139) 


(32) 


Three identification runs were made: 


(a) The underlined parameters in equations (27) and (28) are estimated. 
Note that H is a linear combination of the first two rows of F. The 
estimated and the measured time histories are shown in figure 7. 

(b) The input/output relation is represented as a five pole, four zero 
transfer function. The coefficients of the transfer function are 
identified directly. This requires estimation of 10 parameters. The 
measured and estimated time histories of the accelerometer response 
are shown in figure 8 . The match between the two time histories is 
comparable to that in figure 7. The identified transfer function is 


-0. 0137 


(S^ + 0.432s^ + 0.593s^ + 0.229s + 0.00844) 

(s^ + 0.525s'^ + 0.451s^+ 0.145s^ + 0. 029 7s + 0.00116) 


(33) 
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Certain identifiability problems were indicated by the analysis of 
the information matrix, This^ together with the fact that one zero 
is very close to a pole (eq* (30) ), indicated that a lower order 
model may be more useful, 

(c) The input/output relationship is represented by a four pole, three 
zero model. The coefficients of the fourth order transfer function 
are identified by using the maximum likelihood method. The time 
history plots are shown in figure 9, and the identified transfer 
function is as follows: 


-0,0157 



(s^ + 0.250s^ + 0.569s + 0.139) 

+ 0.467s^ + 0.400s^ + 0.119s + 0.0195) 


(34) 


Though the fit to the time history responses is poorer than before, 
no identifiability problems are indicated, Implying that the fourth 
order transfer function is an adequate representation for the accel- 
erometer output/shaker input relationships. This is also indicated 
by the plots of the gains and the phases of the fifth order and 
fourth order identified transfer functions in figures 10 and 11. 


The poles and zeros of the transfer function used in the simulation and 
the identified transfer functions in each of the three cases given above are 
shown in figures 12 and 13. The closeness of poles and zeros in every case in- 
dicates that the poles and zeros of the transfer function can be identified 
quite accurately and that often lower order models may give as good or better 
results. It shoiild also be noted that the mode which is poorly damped is iden- 
tified more accurately than the mode which has higher damping. 


APPLICATION TO LARGE AEROELASTIC STRUCTURES 


A major application of the input design and instrimient selection procedures 
is the estimation of the flutter characteristics of large three dimensional 
structures. In general, the dynamic flutter behavior of such structures is 
described by partial differential equations in space and time. For the evalua- 
tion of dynamic stability and structural loads, however, a modal analysis or a 
finite element analysis is sufficiently accurate and converts the more complex 
partial differential equations into ordinary differential equations. The dif- 
ferential equations describing the flutter characteristics of large structures 
can be written as follows: 


Mx + Cx + Kx = Du 


(35) 
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where M is the generalized mass matrix, C is the damping matrix, K is the spring 
constant matrix, and D is the input or force distribution matrix* In the 
vector X, each term corresponds to a linear or angular displacement of one of 
the modes. The vector u is a vector of deterministic or random inputs. There 
may be unknown parameters in M, C, K, and D. An accelerometer measures acceler- 
ation at one point, which, in general, is a linear combination of several com- 
ponents of X, i.e., 

(36) 


where represents the parameters which determine the locations of the accel- 
erometers. The strain gages measure displacements, so that their outputs are 
linear combinations of x, as follows: 

Again j determines the locations of various strain gages. The locations of 

the inputs are determined by certain parameters in the matrix D. The transfer 

function between y and u is 
•^a 

y (s) = s^H (Ms^ + Cs + K)"^D u(s) (38) 

a a 

and the transfer function between y and u is 

y (s) = H (Ms^ + Cs + K)“^D u(s) (39) 

s s 

This multivariable problem can be solved in much the same way as described and 
demonstrated above. It is currently being applied to a 19-mode model of tilt 
rotor vehicle for tunnel test. 


CONCLUSIONS 


This paper describes how techniques of system identification can be applied 
to the problem of determining accurate estimates of certain critical parameters 
in aeroelastic structures. 

It was shown that information theoretic approaches can be used for the 
selection of input spectra and instruments. In particular, various functionals 
of the inverse of the information matrix provide useful measures of the accuracy 
of the parameter estimates. The parameter estimation accuracy depends upon both 
the input spectra and the points where the inputs are applied. The functionals 
of the inverse of the information matrix can be minimized under practical con- 
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straints of total costj the cost versus the accuracy of various available in- 
struments and certain location constraints^ the set of instruments which must 
be usedj and the positions at which they should be located can be optimized* 

In some cases » the instruments are already installed^ but the maximum number of 
channels is fixed because of telemetry^ recordings computer capacity^ and many 
other reasons* A tradeoff can be made between various instruments and the 
possibility of having two or more instruments share one channel* 

In the application of system identification methodology to aeroelastic 
structures, the importance of an efficient parameter identification program 
cannot be overestimated. A good method based on the maximum likelihood 
approach that utilizes an eigenvalue analysis of the information matrix provides 
not only efficient parameter estimates but also important diagnostics regarding 
the identifiability of various model parameters and the relevance of various 
models. An aeroelastic analysis of a two dimensional wing, for example, showed 
that a lower order transfer function may be adeqxiate to represent the input/ 
output relationships. 

Maximum utilization of these techniques will, of course, be realized in 
large aerodynamic structures where identifiability could be a serious problem. 
Further work in the development of these techniques is under way. 
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APPENDIX A 


ALGORITHM FOR CHOOSING OPTIMAL INPUT SPECTRUM 


The following procedure may be used to determine the optimal input spec- 


trum; 


(a) Choose a nondegenerate input (i.e., consisting of more than 

frequencies, with a finite power in each frequency). 

(b) Compute the function i|^(a),f) where 


ip(u},f) = Re[Tr R ^ M ^(f) ] to minimize |m ^(f) 


= Re[Tr R 


-1 3T -2 


3T*, 


.- 1 . 


gg WM (f) ] to minimize Tr(M (f)W) 


(Al) 

(A2) 


Find its maximum at w under the constraint of equation (13) . 

o 


(c) Evaluate the information matrix M(w^) at 


(d) Update the design 


f, = (1-a )f+af((o) 0<a<l (A3) 

1 o o o o o • 

is chosen to minimize ]m ^(f) | or Tr(M ^(f)W), where 

M(f-) = (1-a ) M(f ) + a M(w ) , 0 < a <1 (A4) 

1 o o o o o 


It can be shown that such an a exists* 

o 

(e) Repeat steps (b) through (d) until desired accuracy is obtained* 

In the procedure described above^ the function ijj has many local maxima* 

It is computationally time consuming to find O)^ where i|;(a)sf) is maximum* In 

the computer implementation of the algorithm^ we consider finite number of val- 
ues 0 )^ and search through all values of to find the maximum* Most stable 

systems of interest are low pass filters* Thus^ in most cases it is possible to 
find a subset of ^ where the search need be carried out* 
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There are several ways to see if the input design is close to optimal. 

These methods use the following criteria^ 

(a) The information matrix does not change substantially from one step 

to the next, or the optimizing function does not improve significantly 
with iterations* 

(b) The value of which optimizes the value of the desired function 

approaches zero. In other words j little power is placed at newly 
chosen frequencies, 

(c) The maximum value of the function is not much higher than the 

maximum value for the optimum design (i.e,^ m if we maximize |m| and 
Tr(M ^(f^)W) if we minimize Tr(M ^W) . 

In our implementation, (a) and (b) are used as the termination criteria 
and (c) is used as a cheek. 


162 



APPENDIX B; SYMBOLS 
(All Units Are Metric) 


a Elastic axis position from leading edge 

a. Gain coefficient of lift response to step change in angle of 
attack 

a^ (i=l,2...) Coefficients of transfer function denominator 
A State weighting matrix 

A^,...A^ Coefficient matrices defined by equation (2) 

Time constant of lift response to step change in angle of attack 
=b^c/2V 

b. (i-1,2...) Coefficients of transfer function numerator 
r 


c 

C 

C. 

1 

d 

D 


e 

^f 

e 

o 

E 

fC^) 

fo(o.) 

F 

¥ 


Wing chord 

Cost function of all Instruments; generalized damping 

Cost coefficient of a single instrument 

Derivative 

Control distribution matrix to measurements; control matrix 
(equation (35) ) 

Position of c.g-, relative to elastic axis 

Distance of force application from elastic axis (positive forward) 

Distance of accelerometer from elastic axis (positive forward) 

Energy constraint 

Nondegenerate input function 

Starting guess of f(o)) for design algorithm 

Intermediate values of f^(oo) 

System dynamics matrix; total section lift, input force 
=F/2TTqS 
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G 

H 

H 

H 

j 


Control distribution matrix 
Measurement distribution matrix 

Accelerometer measurement state distribution matrix 
Strain gage measurement state distribution matrix 

=/=i 


j 

k 

k 

y 

k. 

k 


m 


K 

L 

L 


1 

2 


m 

m 


1 


Cost functional 

Reduced frequency^ k=(A)c/2V 

Vertical spring rate 

Torsional spring rate 

Mass radius of gyration 

Generalized spring rate 

Quarter chord lift 

Three quarter chord lift 

Mass ; number of parameters in 0 

Virtual mass 


M 

M. . 
P 

q 

r 

r 

o 



Information matrix ; generalized mass 
^ X Elements of information matrix 

(ij=l,2, , .) 

Nondimensional distance = 2Vt/c; number of measurement 

2 

Dynamic pressure = 1/2 W ; number of frequencies 

Nondimensional elastic axis position^ r=e/c 

Nondimensional location of accelerometer, r =e /c 

o o 

Nondimensional input force location, r^=e^/c 
=a/c “1/4 
=3/4 -a/c 

(i=l,2»..) Diagonal elements of measurement noise power spectral density 
matrix 
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S(a),0) 
s^(bi,e) 
S2< O),0) 


T(o),e) 


Tj^(a),6) 

T2(a),9) 


Measurement noise power spectral density matrix 
Real part 

Laplace operator ^"^(d/dp) 

Wing area = cb 

State transfer function 

Specification parameter independent part of S(o),6) 
Specification parameter dependent part of S((a),9) 

Time 

Time interval (equation (10)) 

Output transfer function 

Specification parameter independent part of T(a),0) 
Specification parameter dependent part of T(o),0) 

Lift coefficient, L^/2TTqS; control variable 

Measurement noise 

Velocity 

Nondimensional vertical displacement, y/e 
y/c 

Information matrix weighting matrix 
State 

Measurement vector (p x 1) , vertical deflection 

Accelerometer measurement 

Deflection at accelerometer location 

Power at frequency 

Cutoff power 

Interpolation coefficient for update of the input design 
Input locator parameter 
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Y 

0 

9 

X 

/s 

0 

^ (o) , 0 ) 
X 

1 

y 


y . 

1 

P 

4 ) 



0) 


0) 

o 


0) 

y 

0) 

y 

4> 

0).* 

1 


(I) 


Scalar adjustment factor 
Parameter set 

Parameter of the set 
Parameter estimate 
-l+Tr(AS(a),0)S^(a),0)) 

Lagrange multiplier 
Mass parameter 
Lagrange multiplier 
Air density 
Angular displacement 
Section pitch rate, | 

/s 

Standard deviation of 0 . 

1 

Optimization function 
Circular frequency 

Undamped frequency; starting frequency for design algorithm 
Torsional natural frequency = 

Vertical natural frequency - 
=0)yC/2V 
=w.c/2V 

Frequencies within Au) 

Power weighted average of frequencies 



166 



REFERENCES 


1. Baird, E. F. ; and Clark, W. B. : Recent Developments in Flight Flutter Test- 
ing in the United States. Supplement to the Manual on Aeroelasticity, 
Volume 4, Chapter 10. AGARD Rept. 596, 1972. 

2. Rosenbaum, Robert; Survey of Aircraft Subcritical Flight Flutter Testing 
Methods. NASA CR-132479, Aug. 1974. 

3. Waisanen, P. R. ; and Perangelo, H. J. : Real Time Flight Flutter Testing 
Via Z-Transform Analysis Techniques. AIAA Paper 72-784, Aug. 1972. 

4. Cole, Henry A., Jr.; On-Line Failure Detection and Damping Measurement of 
Aerospace Structures by Random Decrement Signatures. NASA CR-2205, Mar. 
1973. 

5. Houbolt, John C.; Subcritical Flutter Testing and System Identification. 
NASA CR-132480, Aug. 1974. 

6. Gupta, Narendra K. ; and Hall, W. Earl, Jr.: Input Design for Identifica- 
tion of Aircraft Stability and Control Derivatives. NASA CR-2493, Feb. 
1975. 

7. Mehra, Raman K. ; Optimal Input Signals for Parameter Estimation in Dynamic 
Systems - Survey and New Results. IEEE Trans, on Automatic Control, vol. 
AC-19, Dec. 1974, pp. 753-768. 

8. Mehra, Raman K. ; and Gupta, Narendra K. ; Status of Input Design for Air- 
craft Parameter Identification. Methods for Aircraft State and Parameter 
Identification. AGARD-CP-172, Nov. 1974. 

9. Gerlach, 0. H. : The Determination of Stability Derivatives and Performance 
Characteristics From Dynamic Manoeuvres. Presented at 38th Meeting of 
AGARD Flight Mechanics Panel, Toulouse, France, May 1971. 


167 




Figure 1. Various steps in aeroelastic testing. 
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Figure 3. Optimal input design for the id 
parameters. 
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Figure 6. Parameter estimation error as a function 
of accelerometer position. 
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4TH ORDER IDENTIFIED TRANSFER FUNGTIGN 



Figure 10. Gain comparisons of fifth order and 
fourth order identified transfer functions. 
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